A Strong-Coupling Approach to the Magnetization Process of Polymerized Quantum 

Spin Chains 
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Polymerized quantum spin chains {i.e. spin chains with a periodic modulation of the coupling con- 
stants) exhibit plateaux in their magnetization curves when subjected to homogeneous external 
magnetic fields. We argue that the strong-coupling limit yields a simple but general explanation 
for the appearance of plateaux as well as of the associated quantization condition on the magne- 
tization. We then proceed to explicitly compute series for the plateau boundaries of trimerized 
and quadrumerized spin- 1/2 chains. The picture is completed by a discussion how the universality 
classes associated to the transitions at the boundaries of magnetization plateaux arise in many cases 
from a first order strong-coupling effective Hamiltonian. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.45.+j, 75.60.Ej 



Quantum spin systems at low (or zero) temperatures 
can exhibit plateaux in their magnetization curves when 
subjected to strong external fields. Such phenomena 
in quasi-one-dimensional systems have recently been the 
subject of intense interest. In one dimension, there is 
an intriguing interplay between theoretical progress on a 
systematic understanding of the underlying mechanisms 



(see e.g. | 
(see e.g. | 



) and an increasing number of experiments 
I ) on materials which are believed to be pre- 
dominantly one-dimensional. 

Here we study polymerized spin-S' quantum spin chains 
in a magnetic field. Their Hamiltonian is given by 



H — JxSxSx+i — h Sx , 



(1) 



where we assume periodicity of the coupling constants 
with period p, i.e. 



Jx — Jx+p • 



(2) 



We will mostly concentrate on spin S — 1/2 and the 
antiferromagnetic regime > 0. The zero-temperature 
magnetization process of the S = 1/2 polymerized chains 
(|l|) was studied in Q using finite-size diagonalization 
and a perturbative bosonization analysis around the case 
of equal coupling constants = J (apart from this 
and the dimerized case, only some trimerized |^,^ and 
quadrumerized Q cases seem to have been studied in 
the literature). Here we wish to complete the picture by 
discussing the 'strong-coupling' limit where at least one 
coupling constant is small with respect to the others, i.e. 

Jxo ^ 0. 

As is known e.g. from studies of spin- ladders [|,|, 
the magnetization process is easy to understand if some 
Jxo = 0. In this limit, the chain (||) decouples into clus- 
ters of p spins. These 'strongly coupled' clusters mag- 
netize independently such that at zero temperature the 



magnetization (M) can only take finitely many values. 
For spin S they are subject to the quantization condition 



p5(l- (A/)) e Z, 



(3) 



with a normalization such that the magnetization has 
saturation values (M) = ±1. This quantization condi- 
tion was obtained (for 5 = 1/2) in Q and is a special 
case of a more general condition written down in [ 3| . In 
particular the latter was also motivated by considering 
a limit in which the system decouples into clusters of 
finitely many spins. In fact, this counting argument is 
completely independent of the internal coupling inside 
the cluster of the p spins. The quantization condition (|^) 
is therefore insensitive to details of the model. However, 
not only the transition values of the magnetic field but 
also the question if a possible plateau is realized even 
in this limit depends on the precise coupling inside the 
cluster. For the linear arrangement (^ and antiferromag- 
netic Jx > {x ^ xq), all values of (M) permitted by (0) 



are indeed realized at Jx 



0. 



Clearly, it remains also to be shown that the quantiza- 
tion condition (^ is indeed valid at generic points in the 
parameter space, not only for the special points where 
some Jxg = 0. This can be supported by series expan- 
sions around the decoupling point, an issue to which we 
shall return below. 

A first property which one can derive for the full in- 
teracting spin- 1/2 system is the upper critical field /luc 
at which the transition to a fully polarized ferromagnetic 
state takes place. For antiferromagnetic J^, > it is 
simply given by the vanishing of the gap for the one- 
spinwave dispersion above the ferromagnetic background. 
The value of h^c is therefore given by the maximal eigen- 
value of the following p x p matrix 



1 



( Jp+Ji -Ji 

-Ji J1+J2 -J2 







-J2 





\-e-'''Jp 



-e^'^Jp \ 




— Jp-i 



(4) 



where k arises from a Fourier transform using the peri- 
odicity (||) (see also [D for a detailed analysis of a related 
special case). For antiferromagnetic J' > 0, the lowest 
energy excitations occur at A; = for p even and ed k — tt 
for p odd if we introduce the momentum by a translation 
of p sites (as in (|^)). 

In order not to get lost in too many parameters, we 
restrict to the same subspace that was also considered in 
Q before we proceed further. We will now concentrate on 
the following periodic arrangement of coupling constants 



J. - { 



J' 
J 



for X G pZ, 
otherwise. 



(5) 



Now we return to the computation of the largest eigen- 
value of (|). The case p = 2 is a bit special; the correct 
specialization of (^ to p = 2 reads (with the notation 

(i): 



J' + J 



-J 



-ik 



J' 



-J 
j + j' 



J' 



(6) 



Using (||) for p = 2 and (^) for p = 3 and A dX k — pn 
(modulo 27r) we find 



j + j'. 



uc 4 2 4 
h'ir^^ =-J+\{'^'+ V2 J2 - 2JJ' + J' 



(7) 
(8) 

(9) 



respectively. 

Next we turn to series expansions of the plateau bound- 
aries for p < 4. For the present systems, we expect 
that the sharp steps between the magnetization plateaux 
which are present for J' — 01 J — soften as soon as 
one turns on J, J' > 0, but that nothing further happens. 
This scenario was in fact confirmed by the numerical and 
perturbative analysis around J = J' of Q . 

For p = 2 the only non-trivial plateau is located at 
(M) =0. Its boundary is given by the k = spin gap 
E(p=2)_ Series expansions in J'/ J for th is g ap have al- 
ready been carried out some time ago in up to third 
order and have recently been extended to nineth order in 
[0. Adding a further order to eq. (29) of (in pass- 
ing we have also checked eqs. (28) and (30) loc.cit.) one 
arrives at: 



J 



18997 
1769472 



J' 



11960596181 



4892236185600 
+0 (^") . 



21739 
7077888 



J' 



214359199 
6794772480 



833277779047 



117413668454400 



J 
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Here we have used the abbreviation 

J' 



J 



(10) 



(11) 



in order to make the presentation more compact. 

A few remarks may be in place regarding the method 
used here which is summarized e.g. in section 3 of p2[ . 
Like the method of ijllj it exploits the fact that the lead- 
ing coefficients of the series can be obtained on a finite 
lattice. However, we use recurrence relations for the co- 
efficients and an exact symbolic representation through- 
out the computation while in |]ll| a symbolic result was 
reconstructed from a high-precision numerical computa- 
tion. Presumably, cluster expansion algorithms (see e.g. 
llisl ) would be more efficient than the two aforementioned 
methods, but we prefer a simple-minded approach be- 
cause of the ease with which it can be applied to p > 2 
as well. 

For p = 3, there is a plateau at (M) — 1/3, as one 
infers from the above inspection of the case J' = 0. Its 

lower and upper boundaries {h^c~^^ and h^c2~^\ respec- 
tively) are determined by the k — tt gap of the single-spin 
excitations. Up to fifth order in J', one finds the follow- 
ing series: 



J 



,(P=3) 

'-C2 
J 



8 211 2 



77437 



1312200 



J' 



7606883 4 
188956800 



7188324510751 
269989034112000 
521 



J' 



2 18^ 



6480 



J' 



394169 



2260895171 
"79361856000 
+0 {J^) . 



6998400 
535736196039221 

43198245457920000 



J' 
(12) 



Lastly, for p = 4 all relevant excitations have fc = in 
the antiferromagnetic regime J, J' > 0. For reasons that 
should be obvious from the results we content ourselves 
with second order series for the gap E(p=^^ and the lower 
and upper boundaries of the (M) = 1/2 plateau (/ii^~^' 
and hi^~^\ respectively): 



E(p=4) 



i (1 + V3- V2) - ^ (4 + V6 + V2) J 



1 



(3079V6 - 163960V3 + 28775V2 



132480 

-h276026) _|_ Q (• j3 



32 



384 



12288 



2\/6 + 8V2 + 17 
48 



J 



2 



14837760 (-2385712V6 + 4730320^3 
-4947835 V2 



,(P=4) 



= 1 



1 



7747472) 
2\/2-3^ 21691a/2- 



16 



-J 



31648 



43008 



-J' 

(13) 



With the choice of couphng constants (g) there is a 
second decouphng Hmit, namely J — + 0, which for p > 3 
is not equivalent to the case discussed before. This limit 
is special in that several of the coupling constants (||) 
vanish at the same time. This leads to p — 2 free spins 
in zeroth order in J. These free spins are immediately 
polarized once a magnetic field is applied. Only the two 
spins coupled by J' require a finite magnetic field to po- 
larize. This leads to an (M) = 1 — 2/p plateau whose 
upper boundary is given by /lij^ ~ J' — 0{J). 

At first order in J, one now has to perform degenerate 
perturbation theory for the free spins. It turns out that 
at this order they behave as isolated clusters of p— 2 spins. 
The corresponding transition fields have been tabulated 
in 1^ and are indeed a reasonable first approximation to 
the plateau boundaries for 3 < p < 6 of |4| at large J' . 

It is actually not difficult to obtain expansions in J for 
some plateau boundaries. Poor convergence is, however, 
to be anticipated. In the present case, internal properties 
of the decoupled clusters are computed perturbatively 
(which were already taken care of exactly at zeroth or- 
der in the expansions around J' = 0). This is reflected 
e.g. in the fact that the fundamental excitations start to 
disperse (i.e. depend on k) only in the second order in J. 

At p = 3 we find the following eleventh order series for 
the boundaries of the (M) = 1/3 plateau: 



,(P=3) 



2^ 32^ 256 ^ 256^ 



537329 



-J- 



834121 



24576 32768 
15865989569 in . 

84934656 ^ 



310154551 
7077888 



-J- 



— = 7 H 7 

J ^^2 4^ 



— 7^ H 7^ 

64^ ^ 64^ 



(14) 



1317 

8192' 
3539135 



4199 



196608 



96157 



589824 



J- 



133012373 



28311552" 679477248 ^ ' 

This result is valid irrespective of the sign of J. Indeed, 
we find agreement with the second-order result of Q for 
ferromagnetic coupling J < 0. For J < and J' > 0, 
there is an experimental realization of a trimerized sys- 
tem: 3CuCl2-2dioxane. However, since the coupling con- 
stants of this material are roughly given by J / J' ~ —5, 
the experimental magnetization curve JT^ is far outside 
the range of validity of our series (^4|); actually no mag- 
netization plateau is observed experimentally. 



For p = 4 and antiferromagnetic J, J' > 0, one finds 
the following counterpart of (13): 



p(p=4) 1 



0=4) 



J 



= 1 



1 



1 



64^ 
52929 



4 

201 
256 



J- 



32 



J- 



16384 



180845 
65536 



497 
2048 



J 



11887 



■J- 



8192 

o(^-") 



■J- 



Ap=i) 



1 1 



J 



16 



512^ 

6555 

'65536 



4 

-6 



449 



4096 



16 



62051 
524288 



J 



4096 

-10 



(15) 



For the gap E^^'^^^ we restrict to second order only, since 
the high degeneracy at (M) = for J = starts to 
invalidate our approach at the third order. 

We compare our perturbative results to the L = 24 
numerical data Q in Figs. |^|[ The full lines show 
our results for the upper critical fields /luc, the dotted 
lines denote the series expansions around J' = while 
the dashed-dotted lines in Figs. ||,^ show the expansions 
around J = 0. Crosses denote i = 24 numerical data 
of Q and the diamonds show the magnetic fields asso- 
ciated to the plateau-values of (M) at J' = J . For the 
isotropic chain they are /i = for (Af) = and h = 2J 
for (Af) = 1. The fields associated to (M) = 1/3 and 
(A/) = 1/2 are computed from the Bethe-ansatz solu- 
tion of the Heisenberg chain (see e.g. [||). Since Abelian 
bosonization predicts the plateaux to open for J' 7^ J [Q , 
the diamonds denote the expected ending points of the 
magnetization plateaux. 

The case p = 2 is shown for completeness in Fig. ||. 
Since here the regimes J' < J and J' > J are equivalent, 
we display only the former. Here the dotted line shows 
our tenth order series expansion (^0|) for the excitation 
energy at /c = 0. The overall agreement is excellent, as 
has already been observed for the gap in ||ll[] . 

Fig. ^ shows the next case, p — 3. The series ( p^ are 
in excellent agreement with the L = 24 numerical data 
for J' < J. Even the ending point of the (M) = 1/3 
plateau is reproduced quite well. For J' > J, the upper 
boundary of the plateaux is also reproduced reasonably 
by (|lj), while the agreement for the lower boundary is 
poor despite the length of the series. This is not entirely 
surprising as we have remarked above. In fact, inspec- 
tion of the expression for /ic?"'^'' in (|lj) shows that the 
coefficients get larger with increasing order such that this 
series might actually not converge in the region shown in 
Fig.|. 

This comparison of perturbation theory and finite-size 
diagonahzation is completed with the case p = 4 in Fig. 
^. The expansions (|l^) around J' = compare again 



3 



favourably with the numerical data of |Q| although the 
series are only of second order. Also the series for hi^~'^'' 
in (|l5| ) agrees quite well with the numerical data for J' > 

J, while that for /lii"*"* yields good agreement at least 
at the right boundary of Fig. ^. The small- J series for 
E(p=4) is not even shown, since due to its low order it 
cannot be expected to give sensible results in the region 
of interest. As in the case p = 3, the limited quality of 
the series for J' > J can be expected on general grounds 
and is also indicated by inspection of the value of the 
coefficients in (|l5|). 

The comparison of the series ( |l^ ) and (|l^) with the 
numerical data of is complicated by the fact that the 
latter does not extend into the region of small J ~ only 
the region J' <C J is covered well. We have therefore 
performed some further numerical computations for p = 
3 and p = 4. At sufficiently small J one can then nicely 
verify the series order by order - much in the spirit of 
I prf . In this manner we have verified the lowest five to six 
orders of all series in and (^5|) (a standard numerical 
accuracy is not sensitive to the highest orders). 

It has been pointed out recently by several authors 
(see e.g. [p^-pO[) that the strong-coupling approach can 
be extended to describe the transitions between plateaux 
by an effective Hamiltonian. For the case discussed here, 
one will in general have to retain two states per site in 
first order in J'. These two states correspond to the two 
plateau groundstates at J' = between which we wish 
to describe the transition. If the coupling constants are 
chosen to preserve parity (as is the case e.g. for (^)), sym- 
metry arguments imply that the effective Hamiltonian is 
an XXZ-chain. Hence one can immediately carry over 
some well-known universal properties of the XXZ-chain 
(see e.g. section II of |^] for a review) to polymerized 
spin chains. Firstly, the mapping to the XXZ-chain in 
the strong-coupling limits implies that the exponents of 
the correlation functions at the plateau boundaries are 
given by 

Vz=2, Vxy = ^- (16) 

Secondly, the transitions at the plateau boundaries are 
predicted to be of the DN-PT type i.e. the mag- 

netization as a function of applied field h has a square- 
root behaviour close to the plateau boundaries. The 
same conclusions are obtained by the Abelian bosoniza- 
tion analysis of the limit J' — > J Q. This follows from 
results in the theory of commensurate-incommensurate 
transitions which in addition imply that the expo- 
nents ( |l6|) as well as the DN-PT square-root behaviour 
should be universal. The fact that identical conclusions 
are reached by considering two different limiting cases 
are in agreement with such a universal scenario. 

Before concluding, it should be mentioned that the 
quantization condition (^ may have to be relaxed in 
certain cases. For example, it has been shown (see e.g. 
[^,|^) that an (M) = 1/2 plateau can appear if a next- 
nearest neighbour interaction is added to the dimerized 



chain, (Mwith p = 2 (for generalizations of this situa- 
tion see ||25|] ) . This phenomenon can be also understood 
within the strong-coupling analysis lT5|-|l^ if one goes to 
first order, i.e. beyond the naive decoupling limit J' — 0. 
The crucial role is played by the XXZ anisotropy appear- 
ing in the first-order effective XXZ chain. If this effec- 
tive XXZ anisotropy turns out to be sufhciently large 
(A > 1), a gap opens and translational symmetry is 
spontaneously broken. In this manner one finds a further 
plateau precisely in the middle between the two values of 
(M) predicted by considering just the limit J' = 0. This 
illustrates that p in (||) should be taken as the period of 
the groundstate which in general can be different (i.e. an 
integer multiple) of the period of the Hamiltonian. 

To summarize, we have shown that the study of the 
strong-coupling limit not only provides a simple way 
to understand the magnetization process of polymerized 
spin chains qualitatively, but that also quantitatively 
competitive results can be obtained with moderate effort. 
In some respects, the situation is even nicer than for spin 
ladders The bare series in J' yield good results in 

the entire region J' < J, including the ending-points of 
the plateaux. Such favourable conditions are probably a 
special feature of polymerized spin-1/2 chains, as is the 
fact that here the condition (^ is necessary and suf- 
ficient for the appearance of a plateau at J' ^ J. 

It is straightforward to extend the approach of this 
paper to more general interactions or to the computation 
of other quantities. We are confident that further explicit 
strong-coupling computations will provide a useful tool 
e.g. if new experimental data is to be explained in terms 
of model Hamiltonians. 
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for providing me with their numerical data. I am grate- 
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port with the TMR grant FMRX-CT96-0012. 
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FIG. 1. Gap and transition fields for p = 2. The full line 
shows the upper critical field (j^ . The dotted line is our tenth 
order series expansion ( |lo| ) for the spin gap. Crosses show 
the L = 24 numerical data of Cabra and Grynberg and the 
diamonds denote the magnetic fields h a,t J' — J associated 
to (M) = and (M) = 1 respectively. 
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FIG. 2. Transition fields for p = 3. The full line shows 
the upper critical field (^, the dotted lines the series ( p^ ) 
for small J' and the dashed-dotted lines the series ( p^ ) for 
small J. Crosses show L = 24 numerical data of Cabra and 
Grynberg and the diamonds denote the magnetic fields h at 
J' — J associated to (A/) = 1/3 and (M) = 1 respectively. 
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FIG. 3. Same as F ig, pi but for p = 4. The full line shows 
(^), the dotted lines (p"^and the dashed-dotted lines (|l5|). 
Diamonds denote magnetic fields associated to (M) =0, 1/2 
and 1, respectively. 
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